Ising instability of a Holstein phonon mode in graphene 
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We study the thermal distribution of phonons in a graphene sheet. Due to the two electronic bands 
there are two out-of-plane phonon modes with respect to the two sublattices. One of these modes 
undergoes an Ising transition by spontaneously breaking the sublattice symmetry. We calculate the 
critical point, the renormalization of the phonon frequency and the average lattice distortion. This 
transition might be observable in Raman scattering and in transport properties. 
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Transport measurements have revealed many in- 
teresting properties of the two-dimensional mate- 
rial graphene [iHa] . More recently, Raman scatter- 
ing has become an additional experimental method 
to study physical properties of graphene @-|sj]- It 
provides a deeper insight into the elastic properties 
of the honeycomb lattice of carbon atoms and the 
related electron-phonon interaction. Since the lat- 
ter is important also for some of the unusual trans- 
port properties, these results are of great value for 
a better understanding of graphene. 

Recent experiments on hydrogenated graphene 
(graphane) have revealed additional Raman spec- 
tral lines due to hydrogen atoms Q|, indicating 
that there is an important electron-phonon inter- 
action between the conduction electrons and the 
lattice formed by the hydrogen atoms. These re- 
sults indicate that doping can alter the proper- 
ties of graphene by changing the band structure 
as well as by modifiying the elastic properties due 
to phonons. 

Theoretical studies of electron-phonon interac- 
tion in graphene at zero temperature came to a 
mixed conclusion. Although the electron-phonon 
coupling is remarkably strong in terms of a renor- 
malization of the lattice vibrations, its effect on 
transport properties is rather weak 0, [l^. In 
particular, the out-of-plane optical mode (ZO 
phonon) has almost no effect on transport at all 
11]. Some experiments also provide evidence 
that transport properties are not much affected by 
phonons [121 ] . This requires a better understand- 
ing of the electron-phonon interaction and its con- 
sequences for the transport properties. 

In this paper we will explain that the ZO 
phonon mode can play an important role by spon- 
taneously breaking the chiral (sublattice) symme- 
try of graphene. This effect is accompanied by an 
Ising transition which is associated with a soften- 
ing of the optical phonon mode. 



model: The dynamics of phonons of a honey- 
comb lattice is complex and can be described by a 
number of phonon modes [ljj HH, EH, EH ■ For 
simplicity, we focus here on optical phonons, cou- 
pled to the electrons by a Holstein interaction [ll[ • 
Then in a tight-binding description, electrons cou- 
pled to optical (Einstein) phonons at frequency u>q 
are described by the Holstein Hamiltonian as 
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Here c\ (c r ) are the electron creation (annihila- 
tion) operators and b\ (b T ) are the phonon cre- 
ation (annihilation) operators. Then the unitary 
Lang-Firsov transformation reveals that this 
Hamiltonian describes an effective attraction be- 
tween the fermions and a renormalization of the 
electronic hopping amplitude /i r! r<, reducing the 
hopping rate substantially (polaron effect). The 
attractive interaction competes with the repulsive 
Coulomb interaction, where the latter might be 
dominant. 

The conventional approach to determine the 
properties of a coupled system of phonons and elec- 
trons is based on a self-consistent evaluation of the 
self energy (Migdal approximation) [l6|, [ljj ■ The 
latter provides an effective (or renormalized) en- 
ergy and its imaginary part an effective scattering 
rate. Such a static approximation might be insuf- 
ficient in a two-dimensional system, since it does 
not take into account thermal fluctuations. This 
was already discussed in an experimental study 
of graphene Q. To avoid this problem, we in- 
clude thermal fluctuations in our approach. To 
this end, we replace the phonon operators b r , b\ 
by their quantum average: b r (b r ) = v r and 
b\ « (6J) = v*. In this approximation we can keep 
thermal fluctuations but ignore quantum fluctua- 
tions of the phonons. The electrons, on the other 



2 



hand, are studied in full quantum dynamics. This 
reduces the grand-canonical ensemble at inverse 
temperature /3 to a functional integral with respect 
to thermal fluctuations of the lattice distortions v T 
and a trace with respect to the electrons. It should 
be noticed that only the real part of the phonon 
field v couples to the electrons. After performing 
the trace over the electrons we get the partition 
function 

Tre~ m ps J det(l + e - ph )e- 0So V[v] = Z (2) 
with the phonon dispersion 

So = f5>1- (3) 

r 

It is convenient to introduce a sublattice rep- 
resentation for the tight-binding Hamiltonian h, 
since the graphene unit cell contains two atoms, 
each of them has one 7r-orbital. This gives a 
two-component wavefunction, and the electronic 
Hamiltonian can be expressed by 2 x 2 (Pauli) ma- 
trices as 

h e = _^ \ = h x oi + h 2 a 2 - M^o ■ (4) 

with hi = (i+F)/2, h 2 = i(i - F)/2 and with 
the chemical potential (bare Fermi energy) fi. In 
case of the honeycomb lattice t is a matrix that 
connects nearest-neighbor sites on the lattice 

3 
J'=l 

where r is on sublattice A and r+a.j is on sublattice 
B. 

The electron-phonon coupling reads in this sub- 
lattice representation as (spatial dependence of the 
distortions are implictly assumed) 

/ Va V\ - iv 2 \ 
K -^ = a [y l+ iv 2 v B ) 

= ipao + (t><?3 + a(vi<Ti + v 2 <7 2 ) , 

where ip = a(vA + «b)/2 and <f> = ol{va — u.b)/2. 
The in-plane phonon modes v\, v 2 will be ignored 
in the following because they preserve the chiral 
symmetry of the electronic part of the Hamilto- 
nian. Consequently, they are not able to cause 
symmetry breaking. More interesting are the out- 
of-plane modes ip and 4> because they break the 
chiral symmetry of h e . Then the phonon term of 
Eq. © becomes 



We introduce the rescaled frequency uj = uj^/a 2 . 
Eventually, the combined Hamiltonian reads 

h = h e + h e - P h = h%ai + h 2 a 2 + (ip - fi)a + 4>a 3 . 

(5) 

Z of Eq. ([2]) serves as a generating function that 
allows us to get, for instance, the static electronic 
Green's by differentiation of In Z as 

G r , r , = i|(l + e-^) r >-^2?[0, ^ (6) 

with S = S - fi^ 1 lndet(l + e"^). It is impor- 
tant to notice that e~ ps = e ~ 0So det(l + e"* 3 ' 1 ) is 
a non-negative function. Therefore, e~@ s jZ is a 
probability density for the phonon field, and the 
static one-particle Green's function then can be 
rewritten as an average (—}ph with respect to the 
distribution e~ ps /Z 

G r y={(l + e-P h )-l,) ph . (7) 

Without electron-phonon interaction the distribu- 
tion of lattice vibrations is exp(— /3Sq) such that 
the average lattice distortion vanishes: {4>)ph — 
(*0)p/t = 0. A nonzero (ip)ph presents a shift of the 
Fermi energy, whereas a nonzero {4>) p h would break 
the sublattice symmetry. An interesting question 
is whether or not a coupling to electrons can create 
a nonzero average distortion. From the symmetry 
point of view this should not be the case for cj) 
because the system is invariant under the trans- 
formation — > — <j). This is a consequence of the 
fact that the distribution e~ l3S jZ is invariant un- 
der unitary transformations. Choosing the discrete 
unitary transformation U = (<Ti + a 2 )j\f 7 l we get 

UhW = h\o 2 + h 2 ai + (tp - (i)a - cjxr 3 . (8) 

This means for h in Eq. (J5J) an exchange of hi with 
h 2 and a sign change (j) — > — cj). Since the system 
is assumed to be isotropic, only a change (p — > —cj) 
remains. This is the invariance under a Z 2 (Ising- 
like) transformation. Thus {4>) p h 7^ represents a 
typical problem of spontaneous symmetry break- 
ing. 

In previous calculations only the ip mode was 
taken into account, whereas the <j> mode has been 
neglected by assuming <fr = [HI [ll|. Indeed, 
a solution with cj> = always exists due to the 
symmetry in Eq. ([8]). However, spontaneous sym- 
metry breaking is possible. This would be asso- 
ciated with an instability of the = solution. 
Although the symmetry-breaking solution is not 
available from a simple perturbation theory, the 
latter can be employed to analyze the stability of 
the symmetry-preserving solution. For this pur- 
pose we evaluate the shift of the phonon frequency 
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FIG. 1: Frequency shift S0J33 in units of 1/t vs. inverse 
temperature fit for fi/t = 0,3,6 (from bottom to top) 
from Eq. ©. 



with (j)o = ip and 03 = 0. Now we assume 
that the SP solution tp, is uniform in space and 
time. Thus our SP solution is a mean-field ap- 
proximation. This allows us to diagonalize the 
argument of the trace term in Eq. ()10[) using a 
Fourier representation. Moreover, we approximate 
the Hamiltonian h e by its low-energy Dirac behav- 
ior h e ~ k\<j\ + k^oi now. Then the eigenvalues 
of h are A± = ip — fx ± \/0 2 + i 2 /c 2 and the SP 
equations read 



LUIp 



1 



1 

2^ 



sinh(/3(^i — V')) 



cosh(/3(^ - ip)) + cosh(p3^<p 2 + t 2 k 2 ) 



kdk 



(11) 



WoSji — ¥ ujoSji + Sujji (j = 0: ip mode, j — 3: 
mode) for q = (i.e. at the T point of the Bril- 
louin zone). In second-order perturbation theory 
in a [1, [1(3, El this gives for our model a diagonal 
matrix with Sujqq ~ —const./ (3 and 
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cosh(/3tA) + cosh(/3^) 



1 + cosh(^) 



(9) 



where A = 2-^/7? is the momentum cutoff. Thus 
the renormalization of ujq of the ip mode vanishes 
at low temperatures. The situation is different for 
the mode. Although SUJ33 also decreases with de- 
creasing temperature, it has a non-zero value in the 
limit of zero temperature (cf. Fig. [TJ. Then the 
solution — is stable (unstable) if ujq + 6UJ33 > 
(< 0). An unstable situation with negative renor- 
malized phonon frequency of a different mode was 
also found by Castro Neto et al. [5]. 

In order to find a stable solution for the case of 
luq + (5^33 < we must go beyond perturbation 
theory. One possibility is to extend the second- 
order calculation by a partial resummation of the 
perturbation series to all orders of a. This would 
provide a self-consistent calculation for the phonon 
modes. However, it is more convenient here to use 
a saddle-point (SP) approximation for the 0, ip in- 
tegration in Eqs. ©, ©. This calculation can be 
extended to include fluctuations around the SP so- 
lution. The latter allows us to analyze the stability 
of the SP solution, to identify possible instabilities 
and to determine the phonon dispersion. 

The SP approximation is based on the varia- 
tional equation SS = 0: 

^• = Al n [det(l + e-^)] (j = 0,3) (10) 



vV 2 + t 2 k 2 



sin 



h(/3vV> 2 +t 2 k 2 ) 



kdk 



cosh(^ - ip)) + cosh(/V0 2 + t 2 k 2 ) 2tt 

The integration of the second equation can be per- 
formed and gives 



LU(f> 



1 



2ir(3t 2 



1 111 



cosh(^0 2 + t 2 A 2 ) + cosh(^) 
cosh(/?0) + cosh(/3/x) 



with the renormalized chemical potential fi = fi — 

If the mean-field solution — is unstable we 
have to find a solution 7^ of Eq. ([TBI . This 
can be obtained from a quadratic equation after 
an expansion of Eqs. (|TTj) . (ITBl for small and p,. 
Then the equation for reads 



dS 



with the coefficients 

Q = LO 



= fi0 + T0 3 + C/2 2 = 



(14) 



C = 



2iL 



r = 



■ In cosh (fityfn) 



Ant 2 



1 



with t = tanh(/3<y / 7r). Besides the vanishing 
mean-field solution = there are now the two 
non- vanishing solutions 



1 



<>± = ± — J-Vl-Cp 2 



(15) 
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The power law with exponent 1/2 is the result of 
the mean-field approximation. Due to the Ising- 
like symmetry of Eq. ([5]) , the exponent should be 
1/8 instead of 1/2 [l{|. However, this incorrectness 
may only be important very close to the transition 
point. 

Moreover, in a small vicinity around the Dirac 
point p, = the SP equation (jlip can be linearized. 
For small <f> this leads to a renormalized Fermi en- 
ergy 

M+l/w , . 

^ ~ 1 - ln(2)/7r/3^ 2 ' 1 1 

According to the discussion of Eq. (|9|) , a mean- 
field solution 4> = is unstable when fl + Cp 2 < 0. 
For fit ^> 1 we have 

^-2 1 ln2 PP 2 f-\ i-r\ 

n+c » ~"-wt + ^ + ^- (17) 

It should be noticed here that with decreasing tem- 
perature Q, + Cp 2 decreases at the Dirac point 
p = but it increases sufficiently away from the 
Dirac point. This is similar to the behavior shown 
in Fig. Q] In particular, for u> < 1/^/nt and 
fit ~ oo we have <f> = if p 7^ and <f> 7^ for 
/i = 0. Thus, at sufficiently low temperatures a 
symmetry-broken mean-field solution exists only at 
the Dirac point. Away from the Dirac point there 
is no symmetry breaking at very low temperatures. 

For the case f2 + Cp 2 < we must replace = 
by one of the solutions of Eq. (TTSl) . Inserting the 
proper solution into the expression of the renor- 
malized phonon frequency d 2 S/d(f> 2 , we obtain the 
non-negative value for the <f> phonon mode 

c^ = a 2 0«a 2 |ft + C/2 2 | [l + Q{-VL-Cp 2 )] , 

(18) 

where 0(...) is the Heaviside function. The be- 
havior of the renormalized phonon frequency is de- 
picted in Fig. [21 This frequency vanishes as a func- 
tion of the renormalized Fermi energy at the criti- 
cal point n c = y/—Sl/C, provided that 17 < 0. For 
p 2 < [j, 2 we have spontaneous symmetry breaking 
with (Q)ph 7^ and for p? > fi 2 a symmetric phase 
with ((f)) p h = 0. The softening of the phonons 
at the critical point should be observable in Ra- 
man scattering. Symmetry breaking, on the other 
hand, might be observable in transport measure- 
ments as a metal-insulator transition. However, 



large fluctuations of <fi can destroy the gap even 
for ( < I > )ph 7^ 0, similar to the case of a quenched 
random gap |2fj|422l |. 



The values of the model parameter can be com- 
pared with experimental data and bandstructure 
calculations. The electronic hopping parameter is 
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FIG. 2: Renormalized frequency oj^ of the 4> mode (in 
arbitrary units) vs. Fermi energy fi (in units of Q/C) 
from Eq. ((IB) . 



t « 2.7 eV [5( and the bare phonon frequency is 
wo ~ 0.1 eV [23]. Although there is no reliable 
estimate for a, the parameter u) might be small 
enough to have SI < 0, which is necessary to see 
spontaneous symmetry breaking. Even if this is 
not the case, w can be reduced by doping with 
non-carbon atoms (e.g. with hydrogen, oxygen or 
flour). 

In conclusion, we have found that in the Hol- 
stein model for graphene an out-of-plane phonon 
mode can spontaneously break the sublattice sym- 
metry, leading to a spatially fluctuating gap whose 
mean value is nonzero. This effect is accompanied 
by an Ising-like transition. Although it is not clear 
whether this can be observed in pristine graphene, 
doping with non-carbon atoms will provide condi- 
tions to observe such a transition. 
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